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Abstract 

We investigate a two-dimensional classical A^-vector model with a generic nearest- 
neighbor interaction W{ai-aj) in the large- A^ limit, focusing on the finite-temperature 
transition point at which energy-energy correlations become critical. We show that this 
transition belongs to the Ising universality class. However, the width of the region in 
which Ising behavior is observed scales as along the magnetic direction and as 

1/A^ in the thermal direction; outside a crossover to mean-field behavior occurs. This 
explains why only mean- field behavior is observed for N = oo. 

PACS: 75.10.Hk, 05.50.+q, 64.60.Cn, 64.60.Fr 



1 Introduction 



The two-dimensional Heisenbcrg model has been the object of extensive studies which mainly 
focused on the 0(A'^)-symmetric Hamiltonian 



where «Tj is an A'"-dimensional unit spin and the sum is extended over all lattice nearest 
neighbors. The behavior of this system in two dimensions is well understood. It is disordered 
for all finite j3 [1] and it is described for — cxo by the perturbative renormalization group 
[2-4]. The square-lattice model has been extensively studied numerically [5-10], checking 
the perturbative predictions [11-14] and the nonperturbative constants [15-17]. 

Beside Hamiltonian (1), one can also consider more general interactions. In particular, 
one can consider the most general Hamiltonian with nearest-neighbor couplings given by 



where W{x) is a generic function such that W{2) > W{x) for all < a; < 2, in order to 
guarantee that the system orders ferromagnetically for /? ^ oo. The behavior of model 
(2) for /3 — > oo is similar to that of model (1): the critical behavior is described by the 
perturbative renormalization group, the spin-spin correlation length diverging as /9 — > oo. 

While for (3 ^ oo the behavior is similar, for finite values of f3 the more general Hamil- 
tonian (2) may give rise to a more complex phase diagram. In particular, one may observe 
first-order transitions, as it has been shown rigorously for a large class of functions W{x) in 
Refs. [18, 19]. These rigorous results confirm the conclusions of Refs. [20-23] that found first- 
order transitions in the mixed 0{N)-RP^~^ for N = oo, in a modified RP^~^ for = 3, 
and in systems with W{x) ~ x^ for p large enough for iV = 3 and N = oo respectively. The 
presence of first-order transitions does not violate the Mermin- Wagner theorem [1], since 
the spin-spin correlation length is finite at the transition. However, the model studied in 
Ref. [20] ior N — oo also shows a finite-/? critical point. The specific heat diverges indicating 
that energy fiuctuations become critical. On the other hand, the spin-spin correlation length 
is finite, although nonanalytic, at the transition, so that the Mermin- Wagner theorem [1] 
is not violated. For a long time it has not been clear whether the observed critical point 
was an artifact of the large- iV limit or rather a true transition to be expected also for finite 
values of N [20,24,25]. Recently, Blote, Guo, and Hilhorst [21] studied numerically a class of 
systems with W{x) ~ x^ and N — 3 and found for a specific value of p a critical transition 
with essentially the same features as those observed in Ref. [20] for = oo. This confirms 
the existence of finite-/? critical points where energy-energy correlations are critical, while 
the spin-spin correlation function is still short-ranged. 

What remains to be understood is the nature of the critical point. A simple symmetry 
argument [21] indicates that the transition should belong to the Ising universality class for 
any value of A^, and this is confirmed by the numerical results for = 3. On the other 
hand, the large- A^ analysis [20, 22] finds mean-field exponents. In this paper we wish to 
reconcile these two apparently contradictory results. Indeed, we show that, as soon as large- 
N fiuctuations are considered, the transition becomes an Ising one. The mechanism that 
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works in the large- iV limit is very similar to that observed in medium-range models [26-30], 
with playing the role of the interaction range, li t = {P — l3c{N)) / j3c{N) is the reduced 
temperature and /i is a properly defined scaling field that plays the role of the Ising magnetic 
field, wc find that a universal scaling behavior is obtained by taking t^O, /i^O, A^^oo 
and keeping fixed i = tN, h = hN^^"^. Ising behavior is observed only for t <C 1 and h <^ 1, 
so that the width of the Ising critical region goes to zero as N'^^"^ (except in the thermal 
direction where it goes to zero as A^~^) and no Ising behavior is observed for N — oo. In the 
other limit mean-field behavior is observed instead. 

The analysis of the large-iV corrections is quite complex. As in the standard 1/A^ ex- 
pansion we first derive the effective Hamiltonian for the auxiliary fields. However, at the 
critical point the propagator is singular, forbidding a standard 1/A^ expansion. We single out 
the zero mode and show that the effective Hamiltonian that is obtained by integrating the 
massive modes corresponds to a weakly coupled (f)^ theory. The analysis is somewhat com- 
plicated by the fact that the Z2 symmetry is absent and thus there is no natural definition 
of the thermal and magnetic scaling fields. The correct scaling fields are linear combinations 
of the Hamiltonian parameters that must be computed order by order in 1/A'". A detailed 
field-theoretical analysis of the effective Hamiltonian allows us to compute the 1/N correc- 
tions to the position of the critical point and to express the large- A^" scaling behavior close 
to the transition in terms of Ising crossover functions. 

The paper is organized as follows. In Sec. 2 we review the large- A^ limit for model (2) 
on a square lattice, reporting some basic formulas needed in the paper. Derivations are 
reported in Ref. [22]. In Sec. 3 we determine the linear scahng fields for A'^ = 00 and show 
that they are related by the mean- field equation of state. In Sec. 4 and 5 we determine 
the effective Hamiltonian for the auxiliary fields and for the zero mode. Sec. 6 discusses 
the weakly coupled 0^ theory in the absence of Z2 symmetry, determining the relevant 
renormalization constants and the scaling behavior. Finally, in Sec. 7 we apply these results 
to the large- A'^ limit. We determine the correct scaling fields, the 1/N corrections to the 
position of the critical point and the width of the Ising critical region. In Sec. 8 we explicitly 
give the scaling behavior for (cr^<Ta,_|_^), so that we can compare with numerical results, for 
instance with those of Ref. [21]. Finally, in Sec. 9 we report some numerical results for the 
mixed 0{N)-RP^~^ discussed in Ref. [20] and for the Hamiltonian introduced in Ref. [21]. 
Conclusions are presented in Sec. 10. 



2 The large- limit 

We consider Hamiltonian (2) on a square lattice and normahze W{x) by requiring W'{2) — 1 
so that in the spin-wave limit 

H = ^ J dx d^a ■ d^a. (3) 

The standard A^-vector model corresponds to W{x) = x. In this work we do not specify 
the function W{x); we only assume that it depends on a parameter p and that in the 
plane there is a first-order transition line /3 = Pfip) ending at a critical point (PcPc) that 
will be characterized below. 
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Our discussion strictly applies only to ferromagnetic models, but under standard as- 
sumptions can be extended to RP^~^ models which correspond to functions satisfying 
W{1 + x) — W{1 — x). In this case the Hamiltonian is invariant under the local trans- 



formations (Tx — ZxCTx, Zx = ±1. 



The large- iV limit has been discussed in detail in Ref. [22]. Wc introduce three auxiliary 
fields Xxn, pxn-i and fix in order to linearize the dependence of the Hamiltonian on the spins 
cr and to eliminate the constraint — The partition function becomes 

j Y[[dpx^dXx^] Y[Wxd<Tx] e-^^, (4) 

XjX X 

where 

-4 = - - ^ [Xx^i + Xx^l0^x ■ o'x+n - Xxi^Pxi^ + '2W{pxi^)] + 2 {PxO-l - Px) ■ (5) 

X^l X 

We perform a saddle-point integration by writing 



1 _ 



N 

Px = 'y + ^=Px- (6) 

V iV 



In Ref. [22] we computed the corresponding saddle-point equations. They can be written as 

7 = a{4 + ml)/2, (7) 
a - 2W'{t), (8) 

^ = ^W)-2 + T-4^' 

B,{ml) 

where the parameter mo is related to the spin-spin correlation length = I/tuq and 

where the integral is extended over the first Brillouin zone. 

In Ref. [22] it was shown that generic models may show first-order transitions. This 
happens when, for given /3, there are several values of uIq that solve the gap equation (10). 
Here, we will be interested at the endpoint of the first-order transition line mo = moc, P — Pc- 
At this point we have 

where the derivatives are taken at fixed p. We will only consider the generic case, always 
assuming that the third derivative is nonvanishing at the transition. 
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3 Equation of state and scaling fields 

We wish now to parametrize the singular behavior for /3 — > /3c and p — > Pc- Expanding the 
gap equation (10) near the critical point we obtain 

/3 - /3e = E a,^(p - p,r{ml - m^J-. (13) 

nm 

Because of the definition of /5c we have aoo = 0. Moreover, Eq. (12) imphcs that aoi = ao2 = 0. 
For p — p^we see that ml has the leading behavior 

ml-ml^(^Y\ (14) 



2 2 _ / «10(p-Pc)V^^ .T^v 



while for /3 = /3c, we have 



0^03 

The nonanalytic behavior with exponent 1/3 is observed along any straight line approaching 
the critical point, except that satisfying j3 — Pc — aioip — Pc) = 0. Therefore, the correct 
linear scaling field is 

Uh = (3 - Pc- dioip - Pc), (16) 

and we have ml — mlc ~ uj/^ whenever Uh ^ 0. To include the case Uh = 0, we write the 
general scaling equation 

ml-mlc^Uf/^f{x), (17) 

where f{x) is a scaling function and x is a scaling variable to be determined. Then, we use 
again the gap equation (13). Keeping only the leading terms we obtain 

Uh = a2o{p - Pc)^ + an{p - Pc)u]l^f{x) + aQ^Uhf{xf, (18) 

so that 

a-o'ifi.xf + aii(p - Pc)ul^'^f{x) - a2o(p - Pcful^ -1 = (19) 

Since (p — Pc)'^u^^ — [(p — Pc)u^^^^]'^ul^^^ , the third term can be neglected. Thus, we may 
take (the prefactor has been introduced for later convenience) 

x = an{p-Pc)\uh\-^/^. (20) 

The function f{x) satisfies 

ao3f{x)' + xf{x)-1^0. (21) 

Such an equation is exactly the mean-field equation for the magnetization. Indeed, if we 
consider the mean-field Hamiltonian 

n = -hM + ^M^ + (22) 
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the stationarity condition gives 

-h + tM + -M^ = 0, (23) 

which is solved by M — /i^/^/(t|/i|"^/^), where f{x) satisfies Eq. (21) with aos = u/6. It is 
thus clear that the scaling field (16) corresponds to h, while p corresponds to the temperature. 
Note that this identification is not unique, since only the line if = is uniquely defined by 
the singular behavior. For instance, in the usual Ising case, we could define t' = t + ah 
without changing the scaling equation of state, since t\h\~'^^^ = t'\h\~'^^^ + ah^^^. Since the 
scaling limit is taken with h ^ 0, t ^ at fixed t\h\~^/^, we see that ^ t'\h\~^/^. In 

the Ising case, however, there is exact Z2 symmetry and thus the natural t variable is defined 
so that it is invariant under the symmetry. In our case we could define Ut — p—pc+A{(3 — Pc) 

and fix A by requiring the leading correction on any line (except = 0) to be of order u^^^ 

1/3 

instead of order , recovering in this way an approximate Z2 symmetry. For our purposes 
this is irrelevant and thus we will use p — as thermal scaling field. 

Eq. (21) may have more than one solution. A simple analysis shows that, if aos > 0, 
f{x) is the unique positive solution of Eq. (21) and has the following asymptotic behaviors: 
f(x) ~ \/—x/ao3 for x — > —00 and f(x) ~ 1/x for x — > 00. If aos is negative the same 
results apply to —f{—x). 

Eq. (17) gives the leading behavior. It is also possible to compute the subleading correc- 
tions. We obtain for the leading one 

ml - ml = u'J'fix) + ufg{x) + 0{u,), (24) 

with 

_ a2oa:^ + ai2aiixf{xf + 0-040-11/(3;)"^ , . 

aUx + 3aosf{xy) ' ^''^ 

Finally, let us discuss the singular behavior of the energy. We have 

E = 2W{t,p). (26) 

Such a function is regular in rriQ and p. Since p — Pc ~ x\uh\^/^ and — ml ~ u\[^ ^ the 
leading term is obtained by expanding the previous equation in powers of ml — ml- Thus 



E^2Wc + 2W' 



2 



urf{x) + 0{ur). (27) 

mo=rrioc 



where W = dW{x,p)/dx with x — t and the suffix c indicates that all quantities must be 
computed at the critical point. 



4 The calculation: propagator and effective ver- 
tices 

In order to perform the 1/A^ calculation we integrate out the fields cr, which is straightforward 
since the Hamiltonian is quadratic in these fields. We thus obtain 

J Hidp^^dX,^] H dfx, e-^. (28) 
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where the effective Hamiitonian H, can be expanded in powers of the auxiliary fields. If we 
define a five-component field 



* = (/i, Ai, A2,pi,P2), 



then Ti, can be written as 



(29) 



(30) 



1 1 



n=3 



n 



\ Nn/2- 



-J ■ I MEpO E <U.(Pl'---'P-)*^i(PO---*A„(Pn), 



where the indices run from 1 to 5, 

6?p 



(31) 



a=l 



and the integration is over the first Brillouin zone. The propagator can be explicitly written 
1 



as 



P-'(P) 



1 1^0,0 




— 1^0,1 












2^2,0 




-\13W'^ 











2^0,2 





-\(3W'^ 



















V 





-\I3W''' 







1 



(32) 



where W should always be intended as a function of T(mo), 



COS* g-r COS-' Qy 



(33) 



and = 4(sin^]9a;/2 + sin^ Py/2). 

For p — > 0, by using the algebraic algorithm described in App. A of Ref. [13], it is easy 
to express the integrals Aij{0, ml) in terms of the integrals S„(mQ) defined in Eq. (11) with 
n = 1, 2. Explicitly we have 



^00(0, m^) =^2, 
Ao(0,m2) =yloi(0,m^) 



AuiO, ml) = - J(4 + ml)B, + + 8ml + ^0)^2, 
o o 

A2o{0,ml)^Ao2{0,ml) = - + B2 - - {4: + ml)Bi. 



(34) 



^It is useful to write the Fourier transform of A^; as A(p) — e ^^^"/^ J2x ^ ^^^^x- This makes all vertices 
and propagators real. 
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Vertices are analogously computed. It is easy to check that the only nonvanishing contribu- 
tions for which some Ai is equal to 4 or 5 are 

Vtlip,, . . . , Pn) = VtliPl, ...,Pn) = -PW^^\f). (35) 

If all indices satisfy < 3, then 

<U.(P1'---'P-)^(Ep^)= (36) 



n / '^(^i+l - - Pi)T2T 2^^i(Pi' 



+ permutations > , 



where 

^i(P,q) = l, ^2(p,q) = -cos(g^+pj2), i?3(p,q) = -cos(?y +P2//2). (37) 

The permutations should made the quantity in braces symmetric under any exchange of 
(Pj, Ai) (the total number of needed terms is {n — l)\/2). As already discussed in Ref. [22], 
at the critical point the inverse propagator at zero momentum has a vanishing eigenvalue. 
Indeed, a straightforward computation gives 

detP-^(O) = i^o.detsi, (38) 

where Xo.det is given by 

i^o.det = - J^l^ [^B.W - (1 - m^(8 + ml)B^)W"] , (39) 

and 

Eq. (38) shows that the determinant vanishes at the critical point — hence there is at least 
one vanishing eigenvalue — since there dfi/dm^ — 0. The corresponding eigenvector can be 
written as 

computed at the critical point. Indeed, 

5 o 

Y,{P-')ab{Q)zb = ^^^i(0' 1' 1' 0' 0)^- (42) 

B=l ^ 

Note that there is always only one zero mode. Indeed, at the critical point, we have from 
Eq. (40) 

Bf — B2 

so that we can write at criticality 



- -^^W'^ (43) 



_ Bt[Bl-{S + ml)mlBl] 1 

^0,det- 32(52-^2) W''' ^ ' 
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We have verified numerically that the prefactor of l/W is always finite and negative, so 
that i^To.det is always nonvanishing. Thus, for p = Pc the determinant detP~-'^(0) vanishes as 
{ml — mlY, as the eigenvalue associated with the zero mode, cf. Eq. (42). Thus, there can 
only be a single eigenvector with zero eigenvalue. 

Because of the zero mode, it is natural to express the fields in terms of a new basis. 
For each mg and p, given the inverse propagator P^^(p), there exists an orthogonal matrix 
U (p; ml,p) such that U^P~^U is diagonal. If vaIp; ^O'P) = ^ai(p; ^o;P) is the eigenvector 
that correspond to the zero eigenvalue for p = at the critical point and QaoIp! "^Q,p) = 
C/A,a+i(p; "^O'P)' a = 1, ■ ■ ■ 4 are the other eigenvectors, we define new fields $a(p) by writing 

^a(p) = J] t/AB(p)*B(p) = ^^a(p)0(p) + <5^«(P)<^«(P)' (45) 

B a 

where $ = {(p^tpa)- Eq- (45) defines the fields $a up to a sign. For definiteness we shall 
assume ^^(p) to be such that, at the critical point, 

B 

We do not specify the sign of since it will not play any role in the following. 

The new field (j) corresponds to the zero mode, while the four fields arc the noncritical 
(massive) modes. The effective Hamiltonian for the fields $ has an expansion analogous to 
that presented in Eq. (30) for ^. The propagator fUs(p) of $ is 

Pab{p) = J2 PcD{p)UcA{p)UnB{p), (47) 

CD 

while the effective vertices are related to the previous ones by 

K4",l..,A„(Pl>--->Pn) = Yl ^ii!...,S„(Pl>--->Pn)t^BiAi(Pi;mo,p)---C/B„A„(Pn;mo,p). 

Bi,...,Bn 

(48) 

Note that Pab{p) is diagonal by definition, i.e., Pab{p) = (^abPaa{p)- Relation (38) implies 

Ai(0) ~ si ~ (p-Pc): {ml - mlf (49) 

close to the critical point. In Appendix A we prove a very general set of identities among 
the effective vertices. In particular, we show that at the critical point we also have 

yni(O,O,O) = 0. (50) 

This result will play an important role in the following. 

5 Effective Hamiltonian for the zero mode 

As we have discussed in the previous Section, the propagator P has a zero mode at the critical 
point. Thus, in a neighbourhood of the critical point the standard 1/N expansion breaks 
down and a much more careful treatment of the zero mode is needed. For this purpose, we 
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are now going to integrate out the massive modes, obtaining an effective Hamiltonian for 
the field 0. More precisely, we define 



The effective Hamiltonian ?iefr[0] has an expansion of the form 



(51) 



(52) 



where vertices and propagators also depend on N and have an expansion of the form 



m=0 



^"'(p) = E]^^"^'(p)' 



m=0 



V'("nPl,---,Pn) = E]^^-nPl,---,Pn) 



m=0 



We report here the explicit expressions that we shall need in the following: 
^0 = ^ / E P' -P)^«f (P)' 



P ah 

1 



(53) 

(54) 
(55) 



E^ni(P'-P'q'-q)^«''(q)- 

ab 

E ^i2(P' -P' o)4;.(o)v;£;(o, -q) PcM) - 



abed 



E (P' q. -P - q)nS (P, q, -P - q) Ac(q) A<i(p + q) 

V^J'^(p,q,r) = V;?)(p,q,r), 

V^)^^^ (P, q, r, s) = (p, q, r, s) - 

E ^i2(P' q' ~P ~ q)^n6 - ^)Pab{p + q) + two permutations, 

ab 

V^^^ (P, q, r, s, t) = (p, q, r, s, t) - 

E ^na(P' q^ -P - q)^im(i", s, t, -r - s - t)P„b(p + q) + 9 permutations 



ah 



(56) 
(57) 

(58) 
+ 



E ^ni(P' q' -p - q)^iS(r' p + q^ s + t)v^A'i(s, t, -s - 1) 



^(3) 



(3), 



X 



.abed 



xPabiP + q)-Pcd(s + t) + 14 permutationsj , 



(59) 
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where a, b, c, d run from 1 to 4 over the massive modes. 

The identities presented in App. A allow us to derive several relations among the effective 
vertices at the critical point. We have for p — Pc and ml — m^^ 



\/J=')(0,0,0) = 0, (61) 



d ~ 



K^'^ (0,0,0) 



n 2 



dm. 



^g|>V>r(0.0,0,0). (62) 



Relation (60) clearly shows that the standard IjN expansion fails close to the critical 
point. Indeed, outside the critical point, Pg~^(p) is nonsingular and for large it is enough to 
expand the interaction Hamiltonian in powers of 1/A^. On the other hand, this not possible 
at the critical point. Since also the three-leg vertex vanishes at zero momentum in this case, 
cf. Eq. (61), the zero- momentum leading term is the quartic one. Since the coupling constant 
is proportional to 1/A^, the model effectively corresponds to a weakly coupled (^^ theory. In 
order to have a stable 0'^ theory, we must also have t^Q*-^-'(0, 0, 0, 0) > 0. For a generic solution 
of the gap equations satisfying Eq. (12), this is not a priori guaranteed. Note, however, that 
if V'q^'^^O, 0, 0, 0) < 0, then, for p = pc, we have -Po~^(0) ^ a{ml — m^'^, with a < 0, as a 
consequence of Eq. (62): the propagator has a negative mass for N — oo. We believe — but 
we have not been able to prove — that such a phenomenon signals the fact that the solution 
we are considering is not the relevant one. We expect the existence of another solution of 
the gap equation (10) with a lower free energy. 

The weakly coupled 0*^ theory shows an interesting crossover limit. If one neglects fluctu- 
ations it corresponds to tune p and m^ so that H and P~^(0) go to zero as A^ — > oo, in such a 
way that ^A^ and P~^(0)A^ remain constant. In this hmit, Ising behavior is observed when 
the two scaling variables go to zero, while mean-field behavior is observed in the opposite 
case. Fluctuations change the simple scaling forms reported above and one must consider 
two additive renormalizations. A complete discussion will be reported below in Sec. 6. 

Finally, we wish to change again the definition of the field so that the effective zero- 
momentum three-leg vertex vanishes for all p and m^ in the limit A" — > oo. For this purpose 
we now define a new field 

«X(P) = 0(P) + M(p), (63) 

where a and k are functions of p and m^ to be fixed. The function k is fixed by requiring that 
the large- A^ zero-momentum three-leg vertex vanishes. Apparently, all V"*^"-* contribute in this 
calculation. However, because of Eq. (61), Vq (0,0,0) vanishes at the critical point. As 
we already mentioned the interesting limit corresponds to considering = m^ — m^^ — > 
and Ap = p — Pc — together with A" — > oo. We will show in Sec. 7 that this hmit should 
be taken keeping fixed A^A^^/^ and ApN, so that vj^^ (0,0,0) is effectively of order N'^^^. 
Therefore, the equation defining k can be written in a compact form as 



03 , dnk 



n-3 



7V5/6 ^ Nn/2 

n>4 
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where ~ ttno+ani/N is the contribution of the n-point vertex at zero momentum. Eq. (64) 
can be rewritten as 



n=l 



7V(n-i)/3 V 7VV6 



k 



= -03 



which shows that k has an expansion of the form 



(65) 



(66) 



The leading constant Uq depends only on the three- and four-leg vertex, the constant ki 
depends also on the five-leg vertex, and so on. The explicit calculation gives 



K^'^ (0,0,0) 



Vi^^ (0,0, 0,0^ 



+ 



iV V^o(^)(0,0,0)X^-^ (0,0, 0,0,0) ^ Q(^-i/2)_ 



^2t7(5) 



2 [Fo(^)(0,0,0,0)]3 



(67) 



By performing this rescaling, the k-leg 0-vertex that scales with N as N^''^/'^ (except for 
^ ^ 3) gives a contribution to the m-leg x vertex (of course m < A;) of order jVi-™'/2+(m.-fc)/3^ 
which is therefore always subleading. Taking this result into account we obtain 

^siVW / ~^ ~ q)x(p)x(q)x(-p - q) 

^ITKr / ^^^nP>q>r,-p-q-r)x(p)x(q)x(r)x(-p-q-r) + ... (68) 



where 



H 



^-ak [Po-^(0) + lp-(0) 
ak^ ,7(5) 



TV 



+ 



ak'^ ~{3) ^ Oik^.~, 



%^-(o,o,o)-^yJ^)(o,o,o,o) 



24A^3/2 



Vi^'iO, 0, 0, 0, 0) + O{N-^/^), 



P-\P) = a^P,-\p) + ^Pr^(p) - ^^o'^(0' P' -P) + ^^i'\0, 0, p, -p) 



Vo(')(0,0,0,p, -p)+0(iV-4/3), 



6Ar3/2 



V('\p, q, r) = a3V^,(^)(p, q, r) - ^V,^'\o, p, q, r) + ©(iV-^/^), 

V iV 



l-(^) (p, q, r, s) = a%^'^ (p, q, r, s) + ©(iV" ^3) , 
where a — ao + aiN'"^/^ + 0{N~^/^) is fixed by requiring that 



P-\P) - p-\0) = K{p) ^ p2, 



(69) 
(70) 



for p — > 0. 
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6 Intermezzo: The critical crossover limit 



In this section we present a general discussion of the critical behavior of a generic weakly 
coupled lattice 0^ theory, as is the effective Hamiltonian (68) for the zero mode. In partic- 
ular, we show the existence of a critical limit, the critical crossover limit, that allows one 
to interpolate between the Ising behavior and the mean-field behavior. The arguments pre- 
sented here closely follow the discussion of Ref. [29] for medium-range models. The main 
difference is the appearance of odd (f)^ interactions that require additional renormalizations. 

6.1 General considerations 

We wish now to discuss the critical behavior of a model with Hamiltonian 



^^'^ (P' -P - q)9'(p)^(q) v^(-p - q) 
^^■1 1 Is ^^^^ ^' ^' ~P ~ ^ ~ ^Wip)f{^)f{^)f{-p ^ q - s) (71) 

on a square lattice. Here the integrations are extended over the first Brillouin zone, v^(p) is 
the Fourier transform of the fundamental field which is normalized so that K{p) ^ for 
p — > 0. By properly normalizing u (which is assumed to be positive) we can also require 
l^(^)(0, 0, 0, 0) = 1. We shaU only consider the case in which the three-leg vertex vanishes at 
zero external momenta, i.e. y'^-* (0,0,0) = 0. This is an important simplifying assumption. 
Indeed, this means that the cubic interaction has the form (p^O(p, which is a renormalization- 
group irrelevant term according to power counting. 

If y*^'^)(p, q, r) = for all momenta, there is an interesting critical limit, the so-called 
critical crossover limit [31]. Indeed, one can show that, by properly defining a function 
rc{u), in the limit u ^ 0, t = r — rdu) 0, and — > at fixed H/u and t/u, the n-point 
susceptibility Xn, i-e. the zero-momentum n-point connected correlation function, has the 
scaling form 

Xn^ut--f:y^^{t/u,H/u). (72) 

The scaling function /^^"""(x, y) is universal, i.e. it does not depend on the explicit form of 
K{p) and of V^W(p, q, r, s) as long as the normalization conditions fixed above arc satisfied. 
The definition of rc{u) has been discussed in Ref. [31] for the continuum model and in 
Ref. [29] for the more complex case of medium-range models. It is of order u and is obtained 
by requiring X2 to scale according to Eq. (72) in the critical crossover limit. If rdu) — 
riu + 0{v?'), then at one loop we have 



1 u 

X2 



1 /•\/(4)(o,o,p,-p) 



2 A K{p)+t 



t 

For i ^ 

i^(4)(o,o,p,-p) ^ /• 1 , , 1 



+ 0(«2). (73) 



/. 



K{p)+t 



f 1 1 

/ ^ h constant = — — In i -|- constant, (74) 

JpP^ + t 47r 
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where = 4sin^(pi/2) + 4sin^(p2/2). Thus, wc obtain 

~ i " t2 
Therefore, if we define 



1 , 1 , ^ 
ri mu In — h constant 

Stt Stt u 



+ 0{u^). (75) 



fc{u) = -^hiu, (76) 

STT 

the susceptibihty X2 scales according to Eq. (72). No other correction is needed when con- 
sidering the higher-order contributions and thus Eq. (76) gives an exact nonperturbative def- 
inition of rc{u). Such a result can be easily understood. In the continuum two-dimensional 
theory there is only one primitively divergent graph, the one-loop tadpole, and therefore 
only a one-loop mass counterterm is needed to make the theory finite. The same holds for 
the lattice model in the critical crossover limit, which is the limit in which the lattice theory 
goes over to the continuum model. 

Now, let us consider the contributions of the three-leg vertex. Since this is a renormalization- 
group irrelevant term, the only changes we expect concern the renormalization of the bare 
parameters. Indeed, we will now show that, if one properly defines functions rc{u) and hc{u), 
then in the limit -u — >• 0, t = r — rc{u) 0, and h = H — hc{u) — >• at fixed h/u and t/u, the 
n-point susceptibility Xn has the scaling form (72) with h replacing H and with the same 
scaling functions of the symmetric case. 

We will first prove this result at two loops and then we will give a general argument 
that apphes to all perturbative orders. Notice that it is enough to consider the case h — 0. 
Indeed, if Eq. (72) is vahd ior h — and any n, then 

Xn{h) = —,Xn+m{h = 0)/."^ ^ ^ i^^t"- ™/j;-(t/«, 0)/.^ = 

- E -^fn\T{t/u, 0){t/u)-{h/ur. (77) 

ml 

m=0 

which proves Eq. (72) for all values of h. 



6.2 Explicit two-loop perturbative calculation 

At tree level, all contributions to Xn that contain three-leg vertices vanish because 

0. Therefore, Xn is identical to Xn in the Z2-symmetric (p'^ theory; in particular x^n+i — 0. 

At one-loop order, graphs contributing to Xn are formed by a single loop made of a three-leg 

vertices and of b four-leg vertices, with a + 2b — n. Each of them contributes a term of the 

form 

!^ / ^3"(P)^4(P) .7o^ 

where Vsip) = V^^\0,p,-p) and Viip) = V^^^^O, 0, p, -p). The leading contribution for 
t ^ is obtained by replacing each quantity with its small-p behavior, i.e. V4(p) 1, 
^3(p) ~ p'^: and ^(p) p^, with = 4sin^(pi/2) + 4sin^(p2/2). Then, we obtain a 
contribution proportional to 



Jp{P^ + tY^^' ^^^^ 
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P+Q p+q 



Figure 1: The four topologies appearing at two loops. Dots indicate parts of the graphs with 
additional external legs. 



Now, for t ^ we have 



(^2\a r 1 for 6 = 

^ ~ < Int for 6=1 (80) 



lAf + tr^-' Ui-^ for 5 > 2. 

Therefore, for i — > 0, u — > at fixed t/u, we have 

^n/2-l for 6 = 

f'Xn/u ~ { \nt for 6 = 1 (81) 

^ foj- J, > 2. 

Thus, all contributions vanish except those with: (i) n = 1, a = 1, 6 = 0; (ii) n — 2, 
a — 2, b — 0; (iii) n even, b — n/2, a — 0. Contributions (iii) are those that appear in the 
standard theory without Lp"^ interaction. Let us now show that contributions (i) and (ii) can 
be eliminated by redefining rc{u) and hc{u). Consider first xi- At one loop we have 

t hcju) 1 f Vajp) ^ 

-Xi = 7n= / -77r-\ — : + 0{^/u). (82) 



For t — > we have 

Vs{p) f Vs{p) 



J^K{p)+t J^K{p) 



+ 0{tliit). (83) 
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Thus, if we define 



tfien txi/u ~ 0{tlnt/y/u, i/u) — > in tfie critical crossover limit. 
Now, let us consider the two-point function. At one loop we have 



r t Tciu) 1 



-X2 ~ - - 

u u u 



(85) 



The first one-loop term is the contribution of the tadpole that has to be considered in the 
pure ip'^ theory and which requires an appropriate subtraction to scale correctly, cf. Eq. (76). 
The second one is due to the three-leg vertex and is finite for i ^ 0. Therefore, if we define 

rc(^^) = r7-H + ^^, (86) 

where r^^™™(u) is given by Eq. (76), we cancel all contributions of the </?^ vertex. 

Now, let us repeat the same discussion at two loops, in order to understand the general 
mechanism. There are four different topologies that must be considered, see Fig. 1. 

Topology (a). 

A contribution to Xn is proportional to 



where 2a + 2b + c + d — n. The leading contribution for t — > is obtained by setting 
—p, (I, —q) PS! F^'*^(0, 0, 0, 0) = 1. Then, the integral factorizes and we can use 
Eq. (80). Ignoring logarithmic terms, we see that the corresponding contribution to t^Xn/u 
scales as 

^ ^|y/%(e^.)/2 _ (88) 

Thus, a nonvanishing contribution is obtained only for c + d — 0. Three-leg vertices can be 
neglected in the critical crossover limit. 
Topology (h). 

A contribution to Xn has the form 



'2+1 r r 

— / /[l^(^)(0,p,q,-p-q)]^ 



^n/2+1 



l^4(p)"^4(q)^^4(p + q)'^V:3(p)"V^3(q)^V^^(p + q)^ . 
{K{p) + tY+<i+^ (ir (q) + tY+-+'' {K{p + q) + t)'=+/+i ' ^ ' 

where 2{a + h + c) + d + e + f = n — 2. The leading infrared contribution is obtained by 
approximating all expressions with their small-p behavior. Therefore, we can write for t ^ 



u^/^+^ f f ifYifYif+qy 



^" Jp h (^2 + t)«+<i+l(g2 + tf+e+l^^q^ + t)c+/+l 



(90) 
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Integrals of this type can be easily evaluated. By writing = {p"^ + t) —t in the numerator, 
we obtain integrals 

- f f 

JpJ<i{p'^ + t)'^{q'^ + ty{p + q +ty 

with m > 0, n > 0, r > 0. Then, we can rescale p — > t^^'^p, q — > i^/^q and extend 
the integration over all M^. This is possible since the corresponding continuum integral is 
finite. As a consequence the integral scales as ^2-m-n-r ^^le corresponding contribution 
to t^Xn/U' scales as 

Q-/\n/2-a-.-c-l^Q-/\i.^e^m_ (92) 

Thus, a nonvanishing contribution is obtained only for d + e + / = 0. Three-leg vertices can 
be neglected in the critical crossover limit. 

Topology (c). 

A contribution to Xn has the form 



^n/2+l 



[ I V^('H0,p,q,-p-q)l^(3)(p,q,-p q) 

v^4(p)"v^4(q)^i^4(p + ^rvz{j?YvM''v^iv + q)^ .qQ^ 
(x(p) + tY+d+^ {K{q) + ty+^+^ {K{p + q) + ty+f+^ ' ^ ^ 

where 2(a + b + c) + d + e + f = n — 1 and wc assume without loss of generality a > b > c. 
We can now repeat the analysis performed for topology (b). We replace each quantity 
with its small-momentum behavior. In particular, we can replace ^''^•'(p, q, — p — q) with 

{p'^ + q'^ + p + q ). Then, we rewrite each contribution in terms of the integrals Imnr, cf. 
Eq. (91). However, in this case it is possible that one (and only one) of the indices vanishes. 
If this the case, Imnr factorizes and we can use the one- loop result (80). A careful analysis 

shows that in all cases the integral scales as t~°'~^~'^ for t — except when 6 = c = 0. In 
this case, if a 7^ the integral scales as In while for a = it scales as log^ t. Thus, 
ignoring logarithms all integrals scale as t~^~^~'^. Therefore, the corresponding contribution 
to f^Xn/u scales as 

These contributions always vanish. 

Topology (d). 

A contribution to Xn has the form 

f /• „,12 V^i(p)"V^4(q)^V;(p + q)-V^3(p)'^^3(q)^^V^^(p + q)^ 



where 2{a + b + c) + d + e + f = n, a > b > c. We repeat the analysis done for topology (b) 
and (c). We find that the integral scales as 

for a = 6 = c = 0; 
-a-b-c+i ) l/t for 6 = c = 0, a > 1; ^gg^ 

for c = 0, 6 = 1, a > 1; 
otherwise. 
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Thus, if 6 7^ and c 7^ 0, ignoring logarithms, the contribution to f^Xn/u scales as 

j "/2 ^„/2-a-6-c+l ^ j ^^^^ ^(d+e+/)/2+l ^g7^| 

These contributions therefore always vanish in the critical crossover limit. If, however, 
h — c— then 

^y\n,.-a-,-c ^ (]^Y\i.^e^m (Qg) 

and thus one may have a finite (or a logarithmically divergent if a = 0) contribution for 
d — e = f = 0. Let us now focus on this last case in which a — n/2 > with n even. Let us 
now show that these contributions are canceled by the one-loop counterterm due to rc{u). 
Indeed, at two loops we obtain a contribution to Xn of the form 

u-'f v»Wp) fo.. = o ,,,, 



' V Jp [K{p) + ^ \ r« for 6 > 1 

where 2a + b = n. Thus, contributions to f'Xn/u scale as = t^/"^ and thus vanish 

unless 6 = 0. Thus, for each even n there are two contributions that should be considered: 
one with topology (d) and one associated with the counterterm rc{u). Taking properly into 
account the combinatorial factors, their sum is given by (of course, we only consider here 
the contribution to rdu) due to the three- leg vertices) 



u 



n/2+1 



y(^)(p,q,-p-q)2 V('\0,q,-qy 



(100) 



For a > 1, the subtracted term improves the infrared behavior, and indeed the integral scales 
as hit and is therefore irrelevant in the critical crossover limit. 



6.3 Higher powers of the fields 

It is interesting to consider also Hamiltonians with higher powers of the field (p. The standard 
scaling argument indicates that these additional terms are irrelevant for the critical behavior, 
but in principle they can contribute to the renormalization constants like he and Vc- For this 
purpose let us suppose that the Hamiltonian contains also terms of the form 

^^^*r = E— T^/ ■■■/ ^iPi + ---Pk)V^'\pi,---,PkMPi)---^{Pk). (101) 

Note the particular dependence on the coupling constant u, which is motivated by the large- 
N calculation and is crucial in the argument reported below. 

Let us consider the contributions of these additional terms. At tree level we have an 
additional contribution to Xn given by 

7,('^-2)/2 

Axn^—^V^-\0,...,0) (102) 



18 



Thus, f^Axn/u ~ u^'^ that vanishes as m — > for A; > 4. At one- loop order there are 
additional contributions of the form 

j (x(p) + i)-S"^y3(p)"^ • • •K+2(p)""+% (103) 

where Vfc(p) = V^^\—'p^ p, 0, .... 0) and ^^(A; — 2)ak = n. Proceeding as in the previous 
Section, using Eq. (80), and noting that 'Ylik>i^k > by hypothesis, we obtain 

f'Axn/u - -u"/2-iti-Sfc>3«fc „^«3+iEfc>B(fe-4K. (104) 

Here, possible logarithmic terms have been neglected. Then, since some with /c > 5 
is nonvanishing by hypothesis, we find that this correction vanishes in the crossover limit. 
Therefore, no contribution survives at one loop. The same is expected at any perturbative 
order. 



6.4 The general argument 

The discussion reported above shows that at two loops one can define a critical crossover 
limit with crossover functions that are identical to those of the symmetric theory. This is 
expected to be a general result since formally the added interaction is irrelevant. This result 
can be understood diagrammatically. 

Consider the continuum theory with Hamiltonian 

n = \ Y.{d,^f + + ^0^n0 + (105) 

Given an /-loop diagram contributing to the zero-momentum n-point irreducible correlation 
function, we can write the corresponding Feynman integral D as 

D^u^^^''^'' j ... I /(pi,...,pO (106) 

where and A^4 are the number of three-leg and four-leg vertices respectively. Using the 
topological relations Ni = + N^ + l — 1 and + 2N4 — 2l + 2 = n, where iVj is the number 
of internal lines, it is easy to see that 

/(Api, . . . , Ap,) = A2(^-^^)/(pi, . . . , p,). (107) 

Thus, there are primitively divergent diagrams for any n: those with A^4 = (and corre- 
spondingly A^3 = n + 21 — 2) are quadratically divergent, while those with N4 — 1 (and 
A^3 = n + 2/ — 4) are logarithmically divergent. Analogously, diagrams contributing to the 
second derivatives of the n-point irreducible correlation functions with respect to the exter- 
nal momenta diverge logarithmically when N4 = 0. Thus, the renormalized Hamiltonian 
contains an infinite number of counterterms and is given by 

+ $^[z„(A)0" + Cn(A)r-in0], (108) 

n 
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where A is a generic cutoff (in our explicit calculation A is the inverse lattice spacing 1/a, 
that never explicitly appears in the calculation since we set a — 1). Now, let us show that 
all counterterms except those computed in the previous Section can be neglected in the 
critical crossover limit. Suppose that we wish to compute Z„{A) at / loops. Keeping into 
account that the divergence may be quadratic or logarithmic, we expect the /-loop divergent 
contribution to Xn to be of the form 

-——{a,A' + tP.MAyt)] + P2[ln(AVt)]), (109) 

where Pi{x) and P2{x) are polynomials and ai a constant. Therefore, the contribution to 
t'^Xnl'^ vanishes unless n/2 + / — 2 < 0. The only two cases satisfying this condition (of 
course n > 1 and / > 1) are I = n = 1, I = 1 and n = 2, which are the cases considered 
before. Let us now consider the contributions to Xn,i, which is the first derivative of the 
n-point connected correlation function with respect to the square of an external momentum 
p computed at zero momentum. Since momenta scale as t^^^, in the critical crossover limit 
we should have Xn,i ~ ut~^~^fn{t/u,H/u). The divergent contributions are logarithmic 
(diagrams with A^4 = 0) and therefore we expect an /-loop contribution of the form 

n/2+l-l 

— (P[ln(A7t)] + finite terms). (110) 

t 

Considering t'^'^^Xn,i/u, we see that this contribution always vanishes in the critical crossover 
limit. Therefore, the renormahzation constants Cn(A) can be neglected. Thus, the only renor- 
malizations needed are those that we have considered. Finally, let us show that correlation 
functions computed in the renormalized theory have the correct scaling behavior. Indeed, 
in the renormalized theory diagrams scale canonically with possible logarithmic corrections. 
Therefore, D defined in Eq. (106) scales as u^'^'^'^^^H^'^^xlogs, so that the contribution of 
D to t'^Xn/iJ' scales as 

^ X i X ^^4+iV3/2^1-iV4 ^ ^iV3/2_ 

u t'' ^ ^ 

Therefore, the only nonvanishing diagrams have A^s = 0, confirming the claim that three- leg 
vertices do not play any role. 

6.5 A unique definition for the renormalization functions hc{u) and 
rc{u) 

In this Section we wish to discuss again the definition of rc{u) and hc{u). It is obvious 
that these functions are not uniquely defined, since one can add a term proportional to u 
without modifying the scaling behavior. We wish now to fix this ambiguity by requiring that 
t — h — corresponds to the critical point. 

It is easy to see that no modifications are needed for hc{u). Indeed, with the choice (84) 
one obtains the correlation functions of the symmetric theory and in this case the critical 
point is uniquely defined by /i = by symmetry. The proper definition of rc{u) requires 
more care, since we must perform a nonperturbative calculation in order to identify the 
critical point. For this purpose we will use the fact that in the critical crossover limit the 
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perturbative expansion in powers of u is equivalent to the perturbative expansion in the 
continuum 0^ theory once a proper mass renormahzation is performed. 

In the continuum theory, ii i = tcont/'^cont is the adimensional reduced temperature 
defined so that i — corresponds to the critical point, we have at one loop, cf. Eq. (2.10) of 
Ref. [29], 

^ ^ ^' ^''\3 + 8nD,]+0{rnnP), (112) 



^contX2,cont 



t 8nt^ 



In 



where D2 is a nonperturbative constant that can be expressed in terms of renormalization- 
group functions, cf. Eq. (2.11) of Ref. [29]. By using the four-loop perturbative results of 
Ref. [32], Ref. [29] obtained the estimate D2 = -0.0524(2). It is not clear whether the 
error can really be trusted, since in two dimensions the rcsummation of the perturbative 
expansion is not well behaved due to nonanalyticities of the renormalization-group functions 
at the fixed point [33,34]; still, the estimate should provide the correct order of magnitude. 

The expansion (112) should be compared with the perturbative expansion of X2 in the 
lattice model. We write rc{u) as 



f ^ U U f Vsipf 



+ Au, 



(113) 



where ^4 is a constant to be determined. From Eq. (85) we obtain 



UX2 



u u 



log 



32u 



Stt A - An 



By comparing this result with Eq. (112) we obtain 



1 , 2567r 

— m 

Stt 3 



3 

Stt 



^4(p) 1 
X(p) p2 



0{u^). 



(114) 



(115) 



If we use definition (113) with A fixed by Eq. (115), the critical point corresponds to f = 0. 



6.6 Critical crossover limit in the large-A^ case 

In Sec. 5 we showed that the effective Hamiltonian has the form (71) with u ~ 1/N. However, 

in the large- case the parameters that can be tuned are = rn^ — rri^^ and Ap = p — p^, 
and moreover all quantities, beside r and H, depend on these two variables. We wish now to 
determine the changes, if any, that appear in the previous treatment due to the dependence 
of y('=)(pi, . . . , pfe) and K{p) on A^ and Ap. As we show in the next Section, Aj„ and 
Ap scale respectively as and 1/N, so that we can assume an additional dependence 

on u^^^, i.e. we consider K{p;u^^'^) and \/*^^)(pi, . . . , p^; -u^/"^). Note that, by definition, 
K{0; u^^^) = for all values of u. Moreover, we assume, as in the case of interest, that 
V^'^\0, 0, 0; u^^^) = for all values of u. Following the calculation presented in Sec. 6.2, it is 
easy to realize that the dependence of the vertices on u is irrelevant except in hc{u). In this 
case, Eq. (82) becomes 

t K{u) 1 f v,{vy/^) , ^, 
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Because of the prefactor 1 / \/u we must take here into account the first correction pro- 
portional to u^l^. Thus, if Vz{^-yl^) ^ V"3,o(p) + 'U^/V3,i(p) and K^^-^l'^) ^ Ko{p) + 
u^/^Ki{p), we obtain 



t 

-Xi = 
u 



hc{u) 



u 



20^ ip 



V^3,0(P) + «'/'V^3,l(p) M^/=^1^3,0(P)^1(P) 



Ko{v)+t 



+ 0(h^/^). (117) 



Now, t can be set to zero without generating infrared divergences, neglecting corrections of 
order u^^^^tlnt ~ u^^'^lnM. It follows 



hc{u) = 



^3.0(P) i/3 ^3,l(p)i^0(p)-^3,0(p)i^l(p) 

^o(p) Ko{pY 



(118) 



Eq. (118) represents the only equation in which the explicit dependence of the vertices on 
should be considered (Ap is proportional to u and thus it can always be set to zero). In 



all other cases, we can simply set A^ = Ap = 0. 



7 Crossover between mean-field and Ising behavior 

In this Section we wish to apply the above-reported results for the critical crossover limit to 
the Hamiltonian (68) of the zero mode. If we identify 

1^(4) (0,0, 0,0), (119) 

the effective Hamiltonian (52) corresponds to the Hamiltonian discussed in Sec. 6, apart 
from a rescaling of the three-leg and four-leg vertices. 

We begin by defining the additive renormalization constants. The mass renormalization 
constant is given by, cf. Eqs. (113) and (115), 

r(N) = ^Inf^V— (3 + 87rL'2) 

' 87r V 2567r/ 87r^ " 



u 



( [y(3)(0,p,-p)]2 1^^(0,0, p,-p 



2yp\^(4)(0,0,0,0)ir(p)2 yW(0,0,0,0)i^(p)^p2/' ^ ^ 

where D2 is a nonperturbative constant defined in Ref. [29] (numerically D2 ~ —0.052). As 
discussed in Sec. 6.6, we can compute all quantities appearing in Eq. (120) at the critical 
point and keep only the leading terms for — > 00. We thus obtain 

rJN) = Jilnf^V— (3 + 87rL'2) (121) 
' 87r V2567r/ 87r^ ' ^ ' 



/ |[V^,^'Ho,P,-p)Po(p)f -V;)^'^(0,0,p,-p)Po(p) + ?Jv;)^'Ho,0,^ 



where all quantities are computed at the critical point. 

Analogously, we should introduce a counterterm for the magnetic field: 
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As discussed in Sec. 6.6, such a quantity should not be simply computed at the critical 
point, but one should also take into account the additional corrections of order N~^/^. Since 
k ~ 0{N^/^) in the critical crossover hmit, we have 



K{N) 



a 



2VN 



Vr\0, P, -P) - Xn°n?i/ o'^("' P' 
l^J'^(0,0,0,0) 



(0,P,-P)P0(P)] . (123) 

The second integral should be computed at the critical point, while the first one and the 
prefactor of the second one should be expanded around the critical point. Since, as we shall 
show below, = ml — m^^ ~ A^~^/^ and p — Pc ^ N^'^, it is enough to compute the first 
correction in A^. In practice, we find that the renormalization terms have the form 

rc{N) = -^(rcolnA^ + rci), 

hc{N) = -l^(h^ + Amhci). (124) 
V N 



Once rc(iV) and hc{N) are computed, we can define the scaling variables Xf ~ t/u and 
~ h/u. Choosing the normalizations appropriately for later convenience, we define 



X, ^ ^[p-i(0)-r,(7V)], 
N 

Xh = -[H-K{N)]. (125) 
a 

Now, the critical crossover limit is obtained by tuning p and ml around the critical point in 
such a way that Xt and Xh are kept constant, i.e. p pcrit, — * ""^ocrit) N ^ oo at fixed 
Xt and Xh- Note that here Pcrit and mgcrit correspond to the position of the critical point as 
a function of N (thus Pcrit ~^ Pc and m^ ^^^^ — > m^^ as iV — > oo) and are obtaind by requiring 
Xt = Xfi = Q (cf. Sec. 6.5). 

In order to compute the relation between p, m^ and Xt, Xh we set A^ = m^ — m^^ and 
Ap = p — Pc and expand Eq. (125) in powers of A„j and Ap. In the following we shall show 
that A^ ~ N~^/^ and Ap ~ N~^, so that the relevant terms are 

Xt — N[bto^p + 6(1 ApA^ + bt2Al^)+di\nN + d2 (126) 

+N'/''id3o + dnAJ. (127) 

The constants are obtained by expanding P~^(0) and H around the critical point and by 
using the expansions (124). All quantities are analytic in A^ and Ap and several terms are 
absent because of identities (60), (61), and (62). In particular, a term proportional to A^ 
is absent in the equation for Xt- This is a consequence of Eq. (62). Indeed, we have 



«^ (0) - Po (0) 2^(^)^0,0,0,0) 
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1 



d 



^oJ V^o^^^O, 0, 0, 0) V^"^o 
+ O(A^,7V-i), 



(128) 



where in the last step we have used Eq. (62). 

We wish now to determine the behavior of and Ap that is fixed by Eqs. (126) and 
(127). We assume 

X X 

(129) 



A™, = 



TV"' 



A _ "PO 



where 5^0 and 5pQ are nonvanishing constants and a and (3 exponents to be determined. If 
P <l,Eq. (126) imphes 



-3a 



0{N 



(130) 



which requires P = 3a. Now, consider Eq. (127). The term iV^/^A^ is of order A^3/2-o ^^^^ 
cannot be made to vanish. Therefore, we must have f3 > 1. Considering again Eq. (127), 
it is easy to see that all terms containing Ap cannot increase as fast as N^^'^. Therefore, 
cancellation of the term dsoN^^^ requires a — 1/3. Consideration of Eq. (126) imphes finally 
P — 1. This analysis can be extended to the subleading corrections, obtaining an expansion 
of the form 



The coefficients are given by 
<^mo — ~ 



N NV^ 



1/3 



+ 0(A^-^/^). 



(131) 



mO 



3htnd. 



'to "30 
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mO 



3&t0C^30 

1 



[btod3i-bhi{dilnN + d2)] + 
{hhAo - bhibt2), 



3btodso 



Xt + 



-J- i^mobt2 + di\nN + d2- Xt) . 



(132) 



We are not able to compute 5m2 and 5pi but we can however compute a relation between 
these two quantities. We obtain 



35mo5m2bh2, 



+ 



bhl bhlSmO 

Correspondingly, by using Eq. (13), we can compute the expansion of Uh'. 

_ UhO Uhl Uh2 



(133) 



(134) 
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where Uhi depends on Xt and Uh2 on Xh- We thus define a new scahng field by requiring that 
no term proportional to XtN~'^/^ is present. For this purpose we set 



Xy, 



where the coefficient Xmix is given by 



Xn 



The new scaling field has an expansion of the form 



^'^ - AT + Ar4/3 



+ 



Uh2 

Ar3/2' 



where 



Uhl 
Uh2 



-^[-003^31 + (0046/18 - 0036/14)^^0] 



-Xh- 



(135) 

(136) 

(137) 

(138) 

(139) 

(140) 



Interestingly enough, in 111^2 all terms proportional to the unknown quantity 5m2 cancel. At 
this point we can easily compute the 1/N expansion of the critical point Pcrit, /^crit- It is 
enough to set Xh — Xt — Q the previous expansions, obtaining 



Pcrit 



+ 0(Ar-7/6)^ 



xt=0 



/^crit — /^c + 



Urn + aw5po\xt=o ^ 0{N~'^/^) 



It follows that 



P - Pcrit 



2 6ti6/ii5^o 



+ 1 



Xt 



3 rfso / htoN 

3ao3^to ~ Oii^ti Xh 



(141) 
(142) 

(143) 
(144) 



where W/i^crit is the value of Uh at the critical point. Therefore, the 1/N corrections modify 
the position of the critical point and change the magnetic scaling field which should now be 
identified with Uh — uh — Uh^cvit (as we already discussed the thermal magnetic field is not 
uniquely defined and we take again p — Pcrit)- Eqs. (143) and (144) also indicate which are 
the correct scaling variables. Ising behavior is observed only if N{p — Pcrit) and N^^'^{uh — 
'W/i,crit) are both small; in the opposite case mean-field behavior is observed. Therefore, as 
N increases the width of the critical region decreases, and no Ising behavior is observed at 
N — 00 exactly. 
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8 Critical behavior of {<Jx • cTx+h) 

In this Section we wish to compute the large- A?" behavior of 

E= {(T^-(T^+^). (145) 

Such a quantity docs not coincide with the energy. However, as far as the critical behav- 
ior is concerned, there should not be any significant difference. In order to perform the 
computation, note that the equations of motion for the field A^,^ give 

{p,^) = l + E. (146) 

Thus, we have 

E^T-1 + -^{p.,) = T-1 + ^Y. U4b{0){^b.). (147) 



TV VN ^ 

Therefore, we need to compute the correlations (0) and {(fia)- For the zero mode we have 
immediately 

= a{x.) - k = -fr'^ixt, X,) - k, (148) 

where /i^™™(a;t, x/i) is the crossover function for the magnetization in the Ising model and 
the constant k diverges as N^^^, modulo some obvious normalizations. 

Let us now consider {ipa) and show that such a correlation vanishes as iV — >• oo in the 
critical crossover limit. Indeed, we can write 

i^a) = -4ffPab{0) /[T4ed(0,p,-p)Ad(p) + (\>m(0,p,-p)-nii(0,0,0))Ai(p)] 

V-/V Jp 

+^Pab{0)Vbii{0,0,0){4>l)- (149) 

V J\ 



Note that in the last term we have replaced 

Ai(p) = {<Pl), (150) 



a necessary step, since the integral diverges at the critical point and therefore should be 
computed in the effective theory for the zero mode. Now, we have 

{<l>l)^k'-2ak{x.)+a'{xl). (151) 

In the critical crossover limit, the two expectation values are replaced by the crossover 
functions for the magnetization and the energy and by a regular term. Therefore, for N ^ oo, 
the leading behavior is {(f)^) = /c^ ~ N^^^. It follows that {(pa) ~ N~^/^. In conclusion we 
can write 



—Ji {xt,Xh)-k 
:Breg + -^U,,{0)-fr'^{xt, xn) + 0{N-'l'), (152) 
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where £^reg is the regular part of E: 



V N 



= To-l + [Ti-A;ot/4i(0)]5^oiV-^/' + C'(iV-2/3), (153) 



where we have written r ~ tq + TiAm and A; ~ koAmVN. 

Equation (152) shows that the singular part of E behaves as the magnetization in the 
Ising model. For Xt <^ 1 and Xh <^ 1 one observes Ising behavior and thus 

E — Ej-cg ~ \xtf^ for Xh = 0, low t phase 

E-E,,^r^\xh\^/^' for XhT^O, (154) 

where /?/ — 1/8 and 5/ = 15. On the other hand, in the opposite limit we have 

E — Ereg ~ \xtf'^'' for Xfi — 0, \xt\ — > oo, low t phase 

E -E^^^r^ \xh\^/^^^ for Ix^Hoo (155) 

with Pmf = 1/2 and 6mf = 3. Note that the limit \xt\ oo and \xh\ — * oo should always 
be taken close to the critical limit. Therefore, \xh\ — oo means that we should consider 
N ^ oo, Uh ^ Uh,crit, P Pcrit in such a way that N^/'^{uh - Uh,crit) ^ oo, i-e. N should 
increase much faster that the rate of approach to the critical point. 

9 Numerical results for selected Hamiltonians 

In this Section we present some numerical results for some selected Hamiltonians. First, as 
in Ref. [21], we consider 

W{x) = ^-{^y. (156) 
Second, we consider the mixed 0(A'")-RP^~^ model with Hamiltonian [20] 

n = -NPv ^(o-, • (T,+^) - ^ ^(o-, • 0-,+^)'. (157) 
This Hamiltonian corresponds to the function 



W{x) ^px + ^{l-p)x'^, (158) 

where we set (3v = (1 +p)/3/2 and Pt = (1 — p)/3/2. This Hamiltonian is ferromagnetic for 
p > -1. Note that for p = -1 we obtain the RP^'^ Hamiltonian [20,24,25,35^1] 

n^-^j2(^x-<T,+,r, (159) 
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^1 




MC 


1.334721915850 


0.9181906464057 




4.537856778637 


-0.9707166650184 




0.1501849439193 


0.8657494320430 


Clin 


0.4359516292302 


-1.1359750388653 




0.5042522341176 


-0.6248171602172 


(2l9 


-1.2793594495686 


0.1122363376128 




-1.3015714087645 


-0.0061080440602 


0,C\A 


67.512019516378 


0.2847340288532 


h+(\ 


-0.04261881236908 


0.21169346791995 




-0.0043954143923 


0.05334711842197 




1.85335235385232 


0.00015229016444 


hhr\ 


-0.0077917231312 


-0.664558775414698 


bhi 


0.085888253027 


-0.094120307536470 


bho 


-0.5785580673083 


0.049596987743194 




— n 991 (SQ'^QQQztm 




dl 


0.205 


0.2678 




0.614 


-0.808514 




0.374 


1.1692 




0.0315969 


0.00933367 



Table 1: Numerical estimates for Hamiltonian (156) (T^i) and (158) {li.2). 

that has the additional gauge invariance cr^ CxCx, = ±1- Under standard assumptions, 
the large- analysis should apply also to this last model: the local gauge invariance should 
not play any role.^ 

In the large- limit [22], the first Hamiltonian has a critical point at /5c ~ 1.335 and 
Pc ~ 4.538. By using the numerical results reported in Table 1 we can compute the first 
corrections to the critical parameters. We obtain 

/3c ^ 1.335 + ^(36.127 + 2.093 In iV), (160) 
Pc ^ 4.538+ (87.92 + 4.80 In AT) . (161) 

Note the presence of a In N jN correction due to the nonanalytic nature of the renormalization 
counterterms. The correction terms are quite large, indicating that the large- results are 
quantitatively predictive only for large values of A^. This is not totally unexpected since [21] 
Pc ~ 20 for N — 3, that is quite far from the large- A^ estimate Pc ~ 4.538. If we substitute 
A^ = 3 in Eq. (161), we obtain pc ~ 35, which shows that the corrections have the correct 
sign and gi\'e at least the correct order of magnitude. 

^The irrelevance of the Z2 symmetry for the large-/? behavior of RP^~^ models have been discussed in 
detail in Ref. [24,40,41]. Thus, in spite of the additional local invariance, RP^~^ models are expected to 
be asymptotically free and to be described by the perturbative renormalization group. 
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For N = oo the second Hamiltonian has a critical point at [20] (3c ~ 0.918 and Pc ~ 
—0.971. Including the first correction we obtain 



In this case the corrections are smaller. However, for < 99.4 they predict < — 1, 
although only slightly. This is of course not possible, since for p < —1 the system is no 
longer ferromagnetic. Thus, we expect the transition to disappear for some N — Nc, with 
Nc ~ 100 (this is of course a very rough estimate). Since Pc ~ —1, we also predict RP^"^ 
models to show a very weak first-order transition for large N . 



In this paper we investigate generic one-parameter families of A*"- vector models that show a 
fine of first-order finite-/? transitions. We focus on the endpoint of the first-order transition 
line where energy-energy correlations become critical, while the spin-spin correlation length 
remains finite, in agreement with Mcrmin- Wagner theorem [1]. We show that, at the critical 
point, the standard expansion breaks down, since the inverse propagator of the auxiliary 
fields has a zero eigenvalue. A careful treatment shows that the zero mode, i.e. the field 
associated with the vanishing eigenvalue, has an effective Hamiltonian that corresponds to 
a weakly coupled one-component 0^ theory. Thus, the phase transition belongs to the Ising 
universality class for any A^, in agreement with the argument of Ref. [21]. In Ref. [22] it was 
shown that for N = oo the transition has mean-field exponents. We reconcile here these two 
results. If Ut and Uh are the linear thermal and magnetic scaling fields, in the critical limit 
a generic long-distance quantity O has a behavior of the form 



Ising behavior. In the opposite limit one observes mean-field criticality. Therefore, the width 
of the Ising critical region goes to zero as A^ — > oo and, even if the transition is an Ising 
one for any A^, only mean- field behavior is observed for A?" = oo. The behavior observed 
at the critical point for A^ — > oo rensembles very closely what is observed in medium-range 
models [26-28,30], with A^ playing the role of the interaction range. Our analysis fully 
confirms the conclusions of Ref. [21]. 

From a more quantitative point of view, we give explicit expressions for the critical values 
/3c and Pc and for the nonhnear scahng fields to order 1/A^. Numerical results are given for 
the Hamiltonian introduced in Ref. [21] and for mixed 0{N)-RP^~^ models [20]. 

Finally, note that, even though the model we consider is two-dimensional, the same 
discussion can be done in generic dimension d. Indeed, for any d the zero mode has a Hamil- 
tonian that corresponds to a weakly coupled 0^ theory. Thus, for d < 4, the phase transition 
is always Ising-like. Some changes, however, should be done in the scaling equations. For 
generic d < 4, Eq. (72) should be replaced by 



Pc = 0.9182 --i: (11.062 - 1.437 In AT), 
Pc = -0.9707+ 4 (2.905 - 1.265 In A^). 



(162) 



(163) 
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(165) 
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Also the expressions for the renormahzation constants Tc and he should be changed: for in- 
stance, in three dimensions we also expect contributions from two-loop graphs, as it happens 
in medium-range models [29]. The main result of this paper, i.e. the fact that the width of 
the Ising critical region goes to zero a,s N ^ oo, holds in any d < 4. However, in generic 
dimension d the natural scaling variables are {p — Pcrit)^^^^*-^"'^-' and {uh — 'U/i,crit)A^^^*-^~'^-'- 
Thus, for d > 2 the Ising critical region shrinks faster as N increases: in three dimensions 
as N~'^ and in the thermal and magnetic directions respectively. 



A Identities among effective vertices 

In this Section we wish to derive a general set of identities for the zero-momentum vertices. 
For this purpose we first rewrite the propagator and the vertices by considering the constants 
a, 7, and r appearing in Eq. (6) as independent variables without assuming the saddle-point 
relations. We define the integrals 

, , f cos* flj, cos-^ , , 



(7-«E^cosg/,)^ 

It is easy to verify that, if a and 7 are replaced by their saddle-point values, a — 2W'{f) 
and 7 = (4 ml)W'{f), then 

1 f COS* Q COS'' Q 

In terms of a, 7, and r the propagator is simply obtained by using Eq. (32) and replacing 
Aij{0,ml)/[W'{T)]'^ with ^^^2(0;, 7). Let us now consider a generic n-leg vertex V^"''. .^^^ at 
zero momentum. It is easy to verify that the only nonvanishing terms with Ai — 4 or Ai — 5 
for some i are (in this Section we do not explicitly write the momentum dependence since 
in all cases we are referring to zero-momentum quantities) 

v}%^V^%^-pW^-\r). (168) 

If all Ai < 3, Eq. (36) gives 

vlt.,A. = li-ir^'^'-^^'in - l)!A..3,n(«,7), (169) 

where ^2 (resp. i^) is the number of indices equal to 2 (resp. to 3). These expressions allow 
us to obtain simple recursion relations for the vertices. 



^l,Ai,...,An Ai,...,A„^ 



T/("+l) I T/(n+l) _ ^ T/(») 

T/{n+l) _ 9 („) 
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where a, 7, r, and P are considered as independent variables. These relations also apply to 

(2) ^ 

the case n — 2, once we identify = P^^. 

Let us consider the projector on the zero mode va, cf. Eq. (45). Keeping into account 
that the symmetry of the matrix at zero momentum implies V2 — Vs and V4 — V5, we obtain 

Close to the critical point we can write, cf. Eq. (41), 

VA^ZA + 0{S,) (172) 

with 

where K is such to have z"^ = 1. Now, let us note that at the saddle point we have (we 
now think of a, r, and 7 as functions of itIq and p) 

= 2VF"^ = -^^z, + 0{s,), (174) 



dml K dml 



where we have used Eqs. (40) and (34). Thus, if we define 



we can rewrite 



E.An+1) _ ^ f d-f d da d dr d\ („) 



d 



- ^75^<..,A„ + 0(.i), (176) 



dm, 







where in the last term a, 7, r, and f3 take their saddle-point values in terms of uIq and p. 
To go further we compute the derivative of za with respect to ttIq. Since 



J2pXb^b = 0{s,), (177) 

we have 

^ ■dml ^ dml 
where — d'^P/d{miy also vanishes at the critical point. 
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Since X^^^^ = 1, dzB/dml belongs to the subspace orthogonal to za- Thus, if P^^ is 
the inverse of P~^ in the massive subspace, we obtain 

1^ = - E P-B^zc + 0{s,) = -1 5: P^sV^'^n^cZn + 0{s,). (179) 

BC BCD 



Finally, let us compute 

' dm, 



A, A 



d 



= C^yCl + r^T.^i?aPal>V^^^^^^ (180) 

where the indices a and 6 run from 1 to 4 over the massive subspace. Since S2 is of order 
p — Pc and niQ — m^^ this relation implies an identity only at the critical point. However, for 
n = 2 we can obtain another relation. Since V^^ = P^^ and Pab^b = 0{s\) we do not 
need Eq. (179) and obtain directly 

Vlii = G^A-/ + 0(si). (181) 

The presence of corrections of order si gives rise to two critical-point identities. Since 
Pii ~ '^1 '^^ obtain 

Viii = 0, 

^111 = C^-^Pu\ (182) 



where all quantities are computed at the critical point. 
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